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Abstract 

The Coulombic systems [ape) and {aae), (appe), (aape) and (Li^^Li^'^e) placed in a magnetic 
field B > 10^^ G are studied. It is demonstrated a theoretical existence of the exotic ion (HeH)'^'^ 
for B > 5 X 10^^ G in parallel configuration (the magnetic field is directed along internuclear axis) 
as optimal as well as its excited states Ivr, 16. As for the exotic ion He2^ it is shown that in spite of 
strong electrostatic repulsion of a— particles this ion can also exist for B > 100 a.u.(= 2.35 x 10^^ G) 
in parallel configuration as optimal in the states lag (ground state), liTu,16g. Upon appearance 
both ions are unstable towards dissociation with He^ in the final state but with very large lifetime. 
However, at i? > 10000 a.u. the ion {HeHf+ becomes stable, while at B > 1000 a.u. the ion 
He^ becomes stable. With a magnetic field growth, both exotic ions become more and more 
tightly bound and compact, their lowest rotational and vibrational energies grow. At the edge of 
applicability of non-relativistic approximation, B ~ 4.414 x 10^'^ G, there are indications that three 
more exotic linear ions [H — He — H)^^ , (He — H — He)^^ and even Li2'^ in parallel configuration 
may also occur. 
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The possibility of formation of new exotic chemical compounds which do not exist in 
: ield-free case is one of the most fascinating features of the physics in strong magnetic fields 
n particular, it was conjectured long ago that linear neutral chains of molecules 
(polymer-type molecules) situated along a magnetic line might exist in astro-physically 
strong magnetic fields. They were called the Ruderman chains. In the year 1999, the 
first theoretical indication for the existence of the exotic one-electron molecular ion H^'^ in 
magnetic fields B > 10^^ G was provided P]. It is a sufficiently long-living state for which 
the optimal configuration is parallel (all protons situated along a magnetic line) decaying 
to H2 + p 2^ . Perhaps, the most astonishing result is an observation of the fast increase 
of the ionization energy of the exotic H^'^ ion with a magnetic field growth. Eventually, for 
magnetic fields 5 > 3 x 10^'^ G the binding energy of the exotic Hf^ ion becomes larger than 
the binding energy of the traditional H2 ion. As a result the exotic ion H^'^ becomes stable 
as well as the traditional ion and the hydrogen atom. Later in the year 2000, it was 
indicated the possible existence of the exotic one-electron molecular ion Hl^ for magnetic 
fields B > 10^'^ G a]- It is worth emphasizing that for all studied one-electron molecular 
systems H2 ^ H^^, the optimal configuration for B > 10^^ G is parallel, when all protons 
were situated along a magnetic line (see, in particular, [6]). Very recently, it was announced 
that two a-particle contained ions, (HeH)^'^ and -^62^ can exist at 5 > 10^^ [3] (for a 
discussion and review the status of one-electron systems, see |^). 

Usually, the investigations of the Coulomb systems in a strong magnetic field were justified 
by the fact there exists a strong magnetic field on a surface of neutron stars or some white 
dwarfs. The recent discovery of absorption features at 0.7 KeV and 1.4 KeV in the X-ray 
spectrum of the isolated neutron star 1E1207.4-5209 by Chandra X-ray observatory 131 and 
its further confirmation by XMM-Newton X-ray observatory (see also |11||) motivated 
more profound and extended studies in a strong magnetic field. In particular, it was assumed 
the hydrogenic atmosphere of 1E1207. 4-5209 with main abundance of the exotic ion H^~^ 
[1^ which may explain the above absorption features. There were alread y fo und two other 
neutron stars whose atmospheres are characterized by absorption features (lJ,llJ], which are 
waiting to be explained. These discoveries pushed us to make a detailed study of traditional 
atomic-molecular systems as well as for a search for new exotic chemical compounds in a 
strong magnetic field. 

The goal of present article is to make a detailed investigation of Coulomb systems made 
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from one electron and several protons and/or a-particles in presence of a strong magnetic 
field. It is natural to expect that due to vanishing rotational energy the optimal configura- 
tion with minimal total energy is the parallel configuration: all massive positively-charged 
particles are situated along a magnetic line. We demonstrate that for sufficiently strong 
magnetic fields which appear on the surfaces of the neutron stars the exotic ions {HeH)'^'^, 
He^ can exist, while the ions {H — He — H)^'^, (He — H — He)^^ and Li^'^ may appear 
to exist as well. The consideration is non-relativistic, based on a variational solution of the 
Schroedinger equation. Hence, the magnetic field strength is restricted by the Schwinger 
limit B = 4.414 x 10^'^ G. Also it is based on the Born-Oppenheimer approximation of zero 
order: the particles of positive charge (protons and a-particles) are assumed to be infinitely 
massive. 

Atomic units are used throughout {h=me=e=l) , albeit energies are expressed in Rydbergs 
(Ry). The magnetic field B is given in a.u. with Bq = 2.35 x 10^ G. 



I. {HeHf+ MOLECULAR ION 

Since long ago the hybrid system {HeHy~^, made out of a-particle, proton and electron, 
(ape), was attempted to explore for field- free case 1^ and for the case of a magnetic field of 
moderate strength B = 1 a.u. In both studies no indication to appearance of a bound 
state was observed. Here we will show that for larger magnetic fields B > 10^^ G the exotic 
hybrid ion {HeH)^~^ can exist in parallel configuration as optimal being a sufficiently long- 
living state. It decays to He~^ +p. For all magnetic fields the binding energy of (HeH)'^'^ is 
slightly smaller than the binding energy of the atomic ion He~^ and their difference decreases 
as a magnetic field grows. In parallel configuration the orbital momentum projection on the 
molecular axis is preserved and eigenstates are characterized by magnetic quantum number 
m, for the ground state m = (Perron theorem). 

The Hamiltonian which describes two infinitely heavy centers a and b of charges Zi and 
Z2 situated along the line forming the angle 6 with the 2;— axis, and electron placed in a 
uniform constant magnetic field directed along the z— axis, B = (0, 0, B) is given by 

n = -A-^-^-^+ ^ + ipA + Ap)+A' , (1) 

(for geometrical setting see Fig. 1). In the case of [HeHY^ the charges Z\ = Z = 2 and 
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a (Zj) B 



FIG. 1: Geometrical setting for the system of two charged centers of the charges Zi and Z2 placed 
on a line which forms the angle 9 with magnetic line (inclined configuration), 0' is mid-point 
between a and h. 

Z2 = 1 correspond to the a— particle and the proton as heavy charged centers. The vector 
potential is given by a certain one-parameter family of vector potentials corresponding to a 
constant magnetic field B = (0, 0, B) 



where ^ is a parameter. The position of the gauge center or gauge origin, where 
A{x, y, z) = 0, is defined hj x = y = 0, with z arbitrary. For simplicity we fix 2; = 0. 
The gauge origin O is chosen to be located somewhere along the line connecting the charged 
centers but not necessarily coinciding with the mid-point O' (see Fig. If ^ = 1/2 we 
get the well-known and widely used symmetric or circular gauge. If ^ = or 1, we get the 
asymmetric or Landau gauge. 

To study the ground state of the Hamiltonian (^J) for the case of the (HeH)"^^ we use the 
variational method with the physically relevant trial function in a way similar to what was 
done for the H2 molecular ion in a strong magnetic field (see [^). A general recipe of the 
choice of the trial function is presented in |2^]. The trial function has a form 



A = B{{^-l)y, ex, 0) , 



(2) 



'Atrial = Ailpi + A2^2 , 



(3) 



where 



^1 



^2 



g-ai^r-i--B[/3i.«x2+/3iy(l-0?/2] _^ jt^^-»2r2-Blf32^ix^+P2y{l-Oy^] 
-a32ri-a4r2-B[feCa.2+/33i/(l-«)?/^] 



(4a) 
(4b) 
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FIG. 2: Total energy curve, Et viz. R, for the {HeH)^~^ molecular ion at i? = 1000 a.u. in parallel 
configuration (9 = 0°). 

here ai...4, Pix,iy, P2x,2y, P3x,3y, A, Ai^2 and ^ are variational parameters. The function ipi sim- 
ulates the incoherent interaction of the electron with charged centers, where the parameter 
A "measures" an asymmetry in the interaction of the electron with a and p. On the other 
side, 11)2 describes the coherent interaction of the electron with a and p. Considering the 
internuclear distance i? as a variational parameter we have in total 10 variational parameters 
(a free normalization of the trial function (jHl) allows us to keep fixed one of the parameters 



The result of calculations shows that the total energy surface Et = Et{B, R,9) for 
B > 1000 a.u. has global minimum at = 0° and a finite internuclear distance R = Reg which 
gives rise a valley when 6 begins to vary. For smaller magnetic fields there exists either no 
minimum or at most some irregularity. Hence if the minima exists the optimal configuration 
for fixed magnetic field B > 1000 a.u. always corresponds to zero inclination, 6 = 0° 
(parallel configuration), see for illustration Fig. El Furthermore, for any fixed magnetic 
field there exists a critical inclination 6cr beyond of which the minimum in the total energy 
curve at fixed inclination disappears. It implies that the system {HeH)'^^ does not exist for 
inclinations larger than the critical inclination. For example, for B = 10000 a.u. the critical 
angle 9cr ~ 8°, which is much smaller than 9cr for H2 

The total energy Et, binding energy Eb and equilibrium distance Reg as well as the 
position and the height of the barrier on the total energy curve of the system (ape) in 




5 



9922 



E.j,(He+) 



9920 



0.4 



0.8 
R (a.u.) 



1.2 



FIG. 3: Total energy curve, Et viz. R, for the (HeH)'^^ molecular ion at B = 10000 a.u. 
for different inclinations; the dash line marks the total energy of the He~^ atomic ion, the solid 
horizontal line shows the lowest vibrational energy of (HeH)^'^. Dotted line corresponds to 2/R. 
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FIG. 4: {HeHf'^ at B = 10000 a.u.: the electronic distribution \^{x,y = 0,z)|V / \'^{x,y, z)\^(ff 
(a) and their contours (b) 

parallel configuration for various magnetic fields are presented in Table HI Binding energy of 
{HeHY^ is always larger than the binding energy of the hydrogen atom. This implies that 
ion {HeH)^^ does not decay to H + a. In general, the binding energy of {HeHy~^ grows 
very fast with the magnetic field increase being smaller than the binding (ionization) energy 
of the He~^ atomic ion 2a]. It continues till B ^ 3x 10^^ G when these two binding energies 



coincide. At larger magnetic fields B >3 x 10^^ G the total energy of {HeHY^ becomes 
lower than the total energy of He^ (see Table I). It implies the following picture. At the 
magnetic fields -B < 3 x 10^^ G the ion {HeH)"^^ is unstable towards a decay to He~^ + p. 
However, at -B > 3 x 10^^ G this decay is forbidden and the exotic molecular ion {HeHy~^ 
becomes stable .27.] . 

In Figs. |21Elthe total energy of {HeH)"^^ viz. internuclear distance R is shown at B = 
1000 a.u. and B = 10000 a. u., respectively, as an illustration. At large internuclear distances 
R, the behavior of the total energy is defined by the interaction of He~^ and p. It can be 
modeled by the repulsive 2/i?-interaction term that is displayed by the dotted line in Fig. El 
which is in good agreement with the results of our variational calculations. 

The equilibrium distance Reg shrinks down quite drastically, reducing almost by a factor 
of three from B = 1000 a.u. to 4.414 x 10^^ G (see TableHl). Thus, the {HeHf+ ion becomes 
more compact as the magnetic field increases. The position of the maximum Rmax of the 
total energy curve Et{R) also reduces with magnetic field increase but not so dramatically 
as for Req (see Table H}. Together with the increase of the height of the barrier (see Table HJ, 
it indicates an increase in lifetime of (HeH)^'^ with magnetic field growth. The lifetime 
is always very large. The electronic distribution is always a single-peak one with the peak 
position shifted towards the charge Z = 2. As an illustration the electronic distribution is 
presented for B = 10000 a.u. in Fig. HI At i? > 5 x 10^^ G the lowest vibrational state of 
{HeHY^ can also exist (see Table E] and for illustration Fig. As for rotational states, 
we calculated the energy of the lowest rotational state following the procedure developed in 
[itI i for the H2 molecular ion. The rotational energy grows with a magnetic field increase 
being systematically larger than vibrational ones unlike the field-free case. For example, at 
B = 10000 a.u. the lowest rotational energy is 6.72 Ry while the lowest vibrational one is 
equal to 0.65 Ry only with the height of the barrier 2.481 Ry (see Table I). For magnetic 
fields B > 5x 10^^ G the well keeps the lowest vibrational state but not the lowest rotational 
state. 

Once a conclusion about parallel configuration as optimal for the ground state is drawn 
the existence of the states with different magnetic quantum numbers (other than m = 0) 
can be explored at 6' = 0°. In order to do it, the following trial function was chosen, 

'^tr^al = p''"' e^^^^O = p''"' 6^'"'^ (^-l + ^2) , (5) 
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TABLE I: Ground state of the (HeH)'^^ ion: total energy Et, binding energy Ef, = B — Et and 
the equilibrium internuclear distance Req (a.u.) in a magnetic field; the binding energy E^-^ of 
the atomic He ion is given Rmax (a.u.) is a position of the maximum of the potential energy 



barrier, the height of the potential energy barrier AE'"*"^ ™™ and the lowest vibrational energy 



B 


Req 


Et (Ry) 


Eb (Ry) 


E^^* (Ry) 


Rmax 




K,, (Ry) 


Trial Function 


1000 a.u. 


0.316 


962.635 


37.365 


40.40 


0.424 


0.061 


0.13 


Eq.JSl 




0.320 


962.694 


37.306 




0.415 


0.045 


0.12 


Eq.(3) [7] 


2000 a.u. 


0.239 


1953.082 


46.918 




0.405 


0.428 


0.25 


Eq.jg 




0.240 


1953.172 


46.828 










Eq.(3) [7] 


10^=^ G 


0.185 


4195.693 


59.626 


61.99 


0.375 


1.197 


0.41 


Eq.JS) 




0.186 


4195.830 


59.489 




0.370 


1.104 


0.41 


Eq.(3) [2] 


8000 a.u. 


0.152 


7927.677 


72.323 


73.84 


0.350 


2.216 


0.59 


Eq.jnj 




0.153 


7927.765 


72.135 










Eq.(3) [2 


10000 a.u. 


0.142 


9922.697 


77.303 


78.43 


0.341 


2.651 


0.65 


Eq.O 




0.143 


9922.906 


77.094 




0.338 


2.481 


0.65 


Eq.(3) [7] 


14000 a.u. 


0.129 


13914.675 


85.325 


85.74 


0.330 


3.43 


0.77 


Eq.O 




0.130 


13914.919 


85.081 










Eq.(3) [7] 


4.414 X 10^3 G 


0.119 


18690.121 


92.858 


92.53 


0.318 


4.22 


0.88 


Eq.JSJ 




0.120 


18690.398 


92.581 




0.315 


3.98 


0.87 


Eq.(3) f7] 



where 

i)i = ^^e-"i^"i-^^^ + Aae""^"^-^^^ , (6a) 

2 2 

^2 = Ase-'^^^^'i-^^''^-*^ + A4e-"5^'^i-"«^2-^*^ , (6b) 

(cf. (ISI)-(4)), here ai,,,6, Pi. .a, ^i...4 are variational parameters, Z = 2. Considering the 
internuclear distance i? as a variational parameter we have in total 14 variational parameters 
(a free normalization of the trial function ^ allows us to keep fixed one of the parameters 
A14). The calculations were carried out for the magnetic quantum number m = 0, —1, —2. 

At m = (ground state) the function (jSJ leads to a certain improvement in total and 
binding energies in comparison with (3)-(4) (see Table H}. Other characteristics are changed 
slightly except for the height of the barrier - it increases in ^ 10% for all studied magnetic 
fields indicating the system {HeHY^ is more stable. This increase in the height of barrier 
appears as a consequence of the fact that the value of the minimal total energy decreases 
while the maximum of the barrier remained almost unchanged in comparison with previous 
calculation [3]. On Fig. El the variational parameters of the trial function (jSJ are shown. It 
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TABLE II: The state Ivr of the [HeH)^'^ in a strong magnetic field. Et is the total energy, 
Ef) = B — Et is the binding energy, where E^'^^ is the total energy of the He'^ atomic ion in 
2p_i state, Req (a.u.) the equilibrium internuclear distance. The binding energy E^^^ of the He'^ 
atomic ion and the binding energy E^ of the Hydrogen atom both in the 2p_i state are shown 
The binding energy of E^^ in Ivr^ state (see «2(|) is given for comparison. 



B Et (Ry) E^ (Ry) E^'^ (Ry) Ef (Ry) E^^^ (Ry) 



10^3 G 


0.263 4212.485 


42.834 


45.65 


17.10 


26.136 


8000 a.u. 


0.210 7947.374 


52.626 


55.13 


20.32 




10000 a.u. 


0.195 9943.497 


56.503 


58.83 


21.56 


34.068 


14000 a.u. 


0.175 13937.213 


62.787 


64.76 


23.53 




4.414 X 10^3 


G 0.160 18714.251 


68.727 


70.31 


25.36 


41.090 



is worth emphasizing the smoothness of their dependence on the magnetic field strength. 

Using the trial function (jSj) it was found that the nearest excited state of {HeH)^^ to 
the ground state is the lowest in energy state at m = — 1, the state Ivr (see Table II). The 
state Itt of (HeH)'^'^ is stable towards a decay {HeHy^{l7r) ^ H{2p_i) + a but it can 
decay {HeH)^^^!^) He^{2p_i) +p. It is worth emphasizing that the binding energies of 
[HeHY^{l'K) and He^{2p_i) are very close. In Table III the results of calculations of the 
{HeHY^ in 15 state (m = —2) are shown. As in the states la, Ivr the ion [HeHY^ in the 
state 15 is stable towards a decay to H{3d-2) + « but it can decay to He~^{3d-2) + P- The 
difference in binding energies of {HeH)'^^{15) and He~^{3d^2) is very small. 

Finally, one can draw a conclusion that the exotic ion {HeHY^ can exist in parallel 
configuration for magnetic fields i? > 5 x 10^^ G as a long-living quasi- stationary state 
with the total energy which is slightly higher than the total energy of the atomic ion He^ 
2^. However, [HeHY^ becomes stable at i? > 3 x 10^^ G and its two lowest energy 
electronic states at m = — 1, —2 can also exist. Their binding energies grow with a magnetic 
field increase. However, the total energy difference of the lowest corresponding states of 
{HeHY^ and He^ at the same magnetic quantum number m reduces as magnetic field 
grows. For given magnetic field their binding energies are reduced with magnetic quantum 
number decrease. 

It is worth commenting that we attempted to study the existence of the bound state of 
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FIG. 5: (HeH)'^^: variational parameters of the trial function (jSJ viz. B for the la state (ground 
state) . 




TABLE III: The state 16 of the (HeH)'^^ in a strong magnetic field. Et is the total energy, 
Eh = B — Et is the binding energy, where E^'^^ is the total energy of the He~^ atomic ion in 
3d-2 state, Reg (a.u.) the equilibrium internuclear distance. The binding energy ^ of the He+ 
atomic ion and the binding energy E'^ of the Hydrogen atom both in the 3d_2 state are shown 
The binding energy of eJ^'^ in 15g state (see |2o|) is given for comparison. 



B Reg Et (Ry) E^ (Ry) E^'^ (Ry) Ef (Ry) E^^^ (Ry) 



10^3 G 


0.324 4219.204 


36.115 


38.93 


14.77 


22.19 


8000 a.u. 


0.251 7955.361 


44.639 


47.34 


17.70 




10000 a.u. 


0.231 9951.968 


48.032 


50.63 


18.82 


29.20 


14000 a.u. 


0.206 13946.450 


53.550 


55.94 


20.60 




4.414 X 10^3 


G 0.187 18724.192 


58.787 


60.91 


22.26 


35.41 



the system made from lithium nuclei, proton and electron in the parallel configuration. The 
Hamiltonian is given by (0) at Zi = 3, Z2 = 1. We used the trial function (3)-(^ at Zi = 3. 
It turns out that at 5 > 7x 10^^ G a total energy curve displays a well-pronounced minimum 
at zero inclination. Nevertheless, no definite conclusion can be drawn, since these values of 
the magnetic field strength are quite far beyond a domain of applicability of non-relativistic 
consideration. 

II. THE MOLECULAR ION iJe^+ 

Can two a-particles be bound by one electron in a strong magnetic field? Although, on 
the first sight the question sounds surprising, it was just recently mentioned [iSj that based 
on semi-qualitative arguments the system (aae) can develop a bound state for sufficiently 
strong magnetic fields. Recently, we announced that this system can aX B > 100 a.u. (= 
2.35 X 10^^ G) 01 . Our present goal is to make a detailed analysis of the existence of the 
bound states of the system (aae). 

The Hamiltonian ^ aX Zi = Z2 = Z = 2 describes two infinitely heavy a— particles 
situated along the line forming the angle 6 with the 2;— axis, and one electron placed in 
a uniform constant magnetic field directed along the 2;— axis As a method to explore the 
problem we use the variational procedure. The ground state trial function is the same as a 
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function which was successful 
field for arbitrary inchnation 



used to explore the H2 molecular ion in a strong magnetic 
[2^. It has a form 



'3 , 



with 



= g-ai^(n+r2)g-B[/3ia=5a;2+/3ij,(l-5V] 



_ ^^-a3Zri-a4,Zr2 _|_ ^-a4Zri-asZr2^ g--B[/33^§x2+/33y(l-Oy2] 



(7) 

(8a) 
(8b) 
(8c) 



(cf. (jU), where ai...4, Pix,iy, P2x,2y, P3x,3y, ^i...3 and ^ are variational parameters. It is worth 
mentioning that the trial function 0-® has provided for H2 the accurate results for B > 
10^ G, which are the most accurate at i? > 10^'^ G so far. Considering the internuclear 
distance i? as a variational parameter we end up with fifteen variational parameters (the 
normalization of the trial function (|7j) allows one to keep fixed one of the ^1,2,3 parameters). 
Each function ■j/'i 2 is a modification of the celebrated Heitler-London or Hund-Mulliken 
functions by multiplication by the corresponding lowest Landau orbital (modified to the case 
of inclined configuration), respectively. The functions ipi^2 describe the coherent (incoherent) 
interaction of the electron with charged centers, respectively. In turn, the function ip^ is a 
modified Guillemin-Zener function. It can be also considered as a non-linear superposition 
of the functions iIji^2- 

We performed accurate variational calculations of the total energy for the system (aae) in 
magnetic fields ranging B ^ 10^° G — 4.414 x 10^^ G for different inclinations. The obtained 
results indicate that for a magnetic field B > 100 a.u. (= 2.35 x 10^^ G) the total energy 
surface Et = Et{0, R) displays clearly a global minimum for a ^ = 0° and a finite value of 
the internuclear distance R = Reg. It indicates to the existence of the exotic molecular ion 
He2^ in parallel configuration as optimal (see Table IV). The equilibrium distance shrinks 
drastically as a magnetic field increases making the system more and more compact. The 
total and binding (ionization) energies grow as a magnetic field increases. It is interesting 
to make a comparison of the binding energy of the -ffcg^ molecular ion with the binding 
energy of the He~^ atomic ion E^^^ . For magnetic fields B < 1000 a.u. the binding energy 
of the atomic ion He^ is larger than one of the molecular ion He^ . Therefore, He^ is 
unstable towards the decay He'^ He^ + a. However, for B > 1000 a.u. the relation is 
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TABLE IV: Ground state of He2^-ion: total energy Et, binding energy Eh and the equilibrium 
internuclear distance Req (a.u.) in a magnetic field; the binding energy E^^* of the atomic He'^ 
ion is given Rmax (a.u.) is a position of the maximum of the potential energy barrier, the 
height of the potential energy barrier the lowest vibrational E'^-^ and rotational E^^^ 



energies in Ry. 



B 




Et (Ry) 


Eb (Ry) 


Ef'^ (Ry) 


Ef^ (Ry) 


Rmax 




E'Lm 


E?ot (Ry) 


100 a.u. 


0.780 


83.484 


16.516 


19.11 


10.291 


1.02 


0.033 


0.026 


0.297 


150 a.u. 


0.640 


130.702 


19.298 














200 a.u. 


0.565 


178.455 


21.545 


24.11 












300 a.u. 


0.480 


274.866 


25.134 


27.56 












1012 G 


0.420 


396.864 


28.668 


30.87 


17.143 


0.90 


1.024 


0.100 


0.810 


1000 a.u. 


0.309 


960.732 


39.268 


40.40 


22.779 


0.82 


2.466 


0.169 


1.41 


10" G 


0.193 


4190.182 


65.137 


61.99 


35.754 


0.70 


7.328 


0.366 


3.718 


10000 a.u. 


0.150 


9913.767 


86.233 


78.43 


45.797 


0.62 


12.25 


0.561 


6.476 


4.414 X 10" G 


0.126 


18677.857 


105.121 


92.53 


54.502 


0.58 


17.19 


0.739 


9.751 



inverted: the binding energy of the atomic ion He~^ is smaller than the molecular ion He^'^. 
The above-mentioned decay becomes forbidden and the exotic ion He^'^ becomes stable (see 
Table IV) . It is quite striking that the binding energy of the -^62^ is about twice larger than 
the binding energy of the H2 molecular ion in the whole region of magnetic fields where 
both ions coexist (see Table IV). 

A potential curve of the total energy of He^^ in parallel configuration (^ = 0°) is always 
characterized by the existence of a barrier (for illustration see Fig. EHZI)- A position of the 
maximum of the barrier Rmax is reduced with a magnetic field increase, but reduction is not 
that sharp as for Req. The height of the barrier ^E'^"-^-^^^ grows very fast with increase of 
a magnetic field, much faster than the binding energy (see Table IV). These two facts can 
be considered as an indication that the lifetime of He^'^ grows quickly with magnetic field 
increase before becoming infinite. 

We make a study of the lowest vibrational state and find that this state exists for B > 



100 a.u. (see Table IV). The energy of the lowest vi 
minimum of the potential curve and written in Ry) is 



^rational state (measured from the 



V 



where /i = 2mpi.oton is the reduced mass of the system of two a-particles (considering 
^neutron ~ ""^proton = 1836.15 mg) and K, is the curvature (in a.u.) of the potential energy 
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FIG. 6: Total energy curve, Et viz. R, for the He2^'^ molecular ion aX B = 1000 a.u. in parallel 
configuration (9 = 0°); the dotted-dash line marks the lowest vibrational energy, the solid line 
marks the total energy of the He~^ atomic ion. 
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FIG. 7: Total energy Et viz. the internuclear distance R for the He2^ ai B = 10000 a.u. for 
parallel configuration, 6 = 0° and for 9 = 10° . Total energy of the atomic ion He~^ is shown by 
solid line, the dashed line marks the lowest vibrational energy. 



curve near the minimum. It seems evident that the higher vibrational states can also exist 
for strong magnetic fields. As for rotational states, we calculated the energy of the lowest 



rotational state following the procedure developed in [17[ for the H2 molecular ion. The 
rotational energy grows with a magnetic field increase being systematically larger than vi- 
brational ones unlike the field-free case. For example, at B = 1000 a.u. the lowest rotational 
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(a) 



(b) 




FIG. 9: Electronic distribution |*(a;,y = 0, J \'^{x,y, z)\'^d^'r (a) and their contours (b) for 
the He2^ molecular ion at B = 10000 a.u. 

energy is 1.41 Ry while the lowest vibrational one is equal to 0.169 Ry only with the height 
of the barrier 2.466 Ry; at -B = 10000 a.u. the lowest rotational energy is 6.02 Ry while 
the lowest vibrational one is equal to 0.561 Ry only with the height of the barrier 12.25 Ry 
(see Table IV). For magnetic fields B > 1000 a.u. the well keeps both lowest rotational and 
vibrational states. 
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The electronic distributions for parallel configuration are illustrated in Figs. IHlini for 
B = 100 a. u. and 1000 a. u., respectively. In a domain of magnetic fields where the -ffcg^ 
ion can exist the electronic distribution has two clear pronounced peaks near the positions 
of the charged centers. With a magnetic field increase the peaks become less and less pro- 
nounced. Near the Schwinger limit B = 4.414 x 10^^ G the electronic distribution becomes 
almost uniform between the positions of charges with exponential-type decay outside. It 
reminds a behavior of the electronic distribution for the if^ molecular ion in parallel config- 
uration in the domain B = — 10^^ G. Average transverse size of the distribution coincides 
approximately to the Larmor radius. 

The total energy surface = Et{B, R,6) for B > 100 a. u. has global minimum at 
6 = 0° and some finite R = Reg. Hence, the optimal configuration for fixed magnetic field 
B > 100 a. u. always corresponds to zero inclination, 6 = 0° (parallel configuration) in 
agreement with a physics intuition (see for illustration Fig. [7j). The study shows that for 
magnetic fields B > 100 a. u. the total energy at equilibrium grows with inclination, while 
the equilibrium distance decreases with inclination. We find that for any fixed magnetic 
field there exists a critical inclination beyond of which the minimum in total energy curve at 
fixed inclination disappears. It implies that the system He^^ does not exist for inclinations 
larger than the critical. This behavior is similar to what was observed for the H2 molecular 
ion [^. For example, for B = 10000 a.u. the critical angle 6cr ~ 28°, which is similar to 6cr 
for H+ 0. 

It can be concluded that the He^'^ ion can exist for B > 100 a.u. in parallel configuration 
as optimal. This raises a natural question about the existence of its excited states. Guided 
by the results for H2 where the lowest-lying excited states were the lowest states but for 
magnetic quantum numbers m = —1, —2 (see perform a study of 1(7^, ^'^u,g, 15g,« 

states. We use the trial functions similar to those which were applied to study the excited 
states of H2 j^l- They have a form of (7) multiplied by a factor pl'^le*"^'^. It turns out 
that positive z-parity states lnu, 16g exist (see Tables V-VI), while the states of the negative 
z-parity 1(t„, liTg, 16u are repulsive. 

It is worth demonstrating the plots of the variational parameters of (7)- (8) viz the 
magnetic field. In general, the behavior is very smooth but it is quite surprising that at 
B ^ 1000 a.u. where the ife2^-ion becomes stable the behavior of almost all parameters is 
visibly changed. 
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FIG. 10: Variational parameters of the trial function ©-(S) as a function of the magnetic field 
strength B for the lag state for the Heg'*' molecular ion. 




TABLE V: He^: the results for the excited state IvTu. Et, Ef, = B — Et are the total and binding 
energies, Req (a.u.) the equilibrium internuclear distance. The binding energy E^'^^ of the atomic 
ion He^ in the 2p_i state is given "^^ . 



H 

±J 






E, (Ry) E^-^ (Ry) E^ (Ry) 


ZUU aj.U. 


0.936 


185.971 


14.029 






300 a.u. 


0.743 


283.474 


16.526 






10^2 G 


0.631 


406.525 


19.007 


21.483 


11.902 


1000 a.u. 


0.448 


973.421 


26.579 


28.749 


16.126 


10^3 G 


0.268 


4209.722 


45.597 


45.645 


26.136 


10000 a.u. 


0.203 


9938.461 


61.539 


58.830 


34.068 


4.414 X 10^3 


G 0.168 18706.914 


76.065 


70.312 


41.090 



TABLE VI: ■ the results for the excited state 15g of the molecular ion in a strong magnetic 

field. Et, Ef, = B — Et are the total and binding energies, Reg (a.u.) the equilibrium internuclear 
distance. The binding energy E^^^ of the atomic He^ ion in the 3(i_2 state is given 2^. 



B Reg Et (Ry) E, (Ry) EI;^^^ (Ry) eI^^ (Ry) 



300 a.u. 


1.00 286.611 


13.389 






1012 G 


0.798 410.089 


15.443 


17.853 


9.861 


1000 a.u. 


0.545 978.234 


21.766 


24.127 


13.488 


1013 G 


0.318 4217.409 


37.910 


38.928 


22.194 


10000 a.u. 


0.238 9948.369 


51.631 


50.633 


29.198 


4.414 X 10i3 


G 0.195 18718.732 


64.247 


60.911 


35.407 



Finally, one can draw a conclusion that the exotic ion He^'^ can exist in parallel config- 
uration for magnetic fields B > 100 a.u. At first, it exists as a long-living quasi-stationary 
state with the total energy which is slightly higher than the total energy of the atomic ion 
He+ |30]. Then Hel+ becomes stable at 5 > 1000 a.u. Making a comparison of the total 
energies of all existing one-electron systems made out of protons and/or a-particles (see e.g. 
^) one can draw a conclusion that at S > 1000 a.u. the exotic molecular ion -ffe^^ has the 
lowest total energy! Two lowest energy electronic states of He^'^ at m = —1, —2 can also 
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exist. For given magnetic field their binding energies are reduced with magnetic quantum 
number decrease. To the contrary to {HeHy^, the total energy difference of the lowest 
corresponding states of He^'^ and He~^ at the same m increases as a magnetic field grows. 

III. THE SYSTEMS {H- He- H)''^+, {He - H - He)^+. 

A. The molecular ion Li^'^ 

Using the variational method as for the He2'^ molecular ion we proceed to a question 
about the existence of the Li^'^ molecular ion in a magnetic field in parallel configuration. 
This system contains two Li nuclei and one electron. It is described by the Hamiltonian 
^ at Z = 3 (see Fig. [T] with Zi = Z2 = 3). We used the same variational trial functions 
©-(jHl) as in the study of Hef^. There was no indication to a minimum in the domain of 
applicability of the non-relativistic approximation B < 4.414 x 10^^ G. Hence, it is unlikely 
the ion Li^'^ can exist in this domain. However, beyond of this domain, at i? > 6 x 10^'^ G 
the total energy curve begins to display a minimum. Although a contribution of relativistic 
corrections for this magnetic field is unknown it seems natural to assume that they do 
not change dramatically the situation. Keeping this in mind, we made calculations for 
B = 7 X 10^'^ G. The result is the following. The total energy displays the well-pronounced 
minimum at i? ~ 0.17a.u. and the minimal energy is equal to 29625.18 Ry. This minimum 
is stable towards small inclinations of the molecular axis. The top of the potential barrier is 
located at i? ~ 0.25 a.u. and the height of the barrier is AE ^ 0.7 Ry. For the same magnetic 
field the total energy of the atomic Li ion is EriLi^^) ~ 29597.14 Ry (see ). Therefore, 
the system (Li^^Li^~^e) is highly unstable towards dissociation Li^'^ Li'^'^ + Li'^^ with 
the large dissociation energy ~ 28 Ry. 

B. The {H -He- Hf+ and {He - H - He)^+ molecular ions 

Let us check the existence of bound states of three-center linear one-electron systems 
{appe) and {aape). If they exist, it would be an indication to the existence of the molecular 
ions {H — He — H)^^ and {He — H — HeY^ , correspondingly. Following the experience 
of (for a discussion see j^) it seems natural to assume that if above systems would 
exist the optimal configuration should be linear and parallel configuration, when all massive 
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B 

FIG. 11: Geometrical setting for a problem of three charged centers situated along a magnetic 
line. The charged centers (marked by bullets) are situated on the z-line at distances R± from the 
central charge which is placed at the origin. 



positively charged centers are situated along a magnetic line. Since two centers are identical 
there can appear two physical configurations depended on where the third, non-identical 
particle is situated: either between two identical particles (symmetric configuration) or as a 
side particle (asymmetric configuration). Geometrical setting is presented in Fig. ^2 Below 
we state that following our calculations a symmetric configuration when the non-identical 
particle is in between the identical particles the bound state for each system can exist, when 
the asymmetric configuration never leads to bound state. 

The variational procedure is taken as a method to explore the problem. The recipe of 
choice of trial functions is based on physical relevance arguments and is described in full 
generality in to where the reader is addressed. The ground state trial function is similar 
to one which was successfully used to explore the H^'^ molecular ion in a strong magnetic 
field jsl [4 1 when the permutation symmetry holds for two centers. The trial function for this 
case is of the form 

'^tr^al = ^1^1 + ^2^2 + A3^P3 , (9) 

with 

^2 = (e""^"^^ + e-"3r-3^g-a4r2g-%Sp2 ^-^Q^^ 
_ ^^-asTi-aara _^ ^-aen-asr^^^-arr^-^Bp'^ (lOc) 
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TABLE VII: Ground states of different one-electron systems at B = 4.414 x 10 G. All energies 
are in (Ry) and Rcq in (a.u.), see text. Total energies and E^"" were calculated variationally, 



System 


Et 


Eb 




E« E^^^ 


{H 


-He- H)'-^+ 


18703.29 


79.69 


0.37 














18750.48 18728.48 18727.75 18690.45 18690.12 18677.86 


{He 


-H- Hef+ 


18712.22 


70.76 


0.34 





where ai...7, (3i,f32,f^3, and v4i...3 are variational parameters. Taking the internuclear dis- 
tances Rj- as variational parameters we end up with 14 variational parameters in total (the 
normalization of the trial function allows us to keep fixed one of the ^1,2,3 parameters). 
The functions ipi describes coherent interaction of the electron with charged centers, when 
ip2 describes incoherent interaction of the electron with identical centers. 

Near the Schwinger limit B = 4.414 x 10^'^ G for both systems {appe) and {aape) in 
symmetric configuration (p — a — p) and {a — p — a) the well-pronounced minimum in total 
energy occurs at = /?_ (see Table EH}- In the same time no indication to appearance 
of the minimum for asymmetric configuration {a ~ p — p) and {a — a — p) is found. For 
both systems the discovered potential well is rather shallow, the potential barrier heights 
are 0.34 Ry and 0.26 Ry, respectively. But as a magnetic field increases it deepens quickly 
for both systems. Both systems are quite compact (here the equilibrium distance is defined 
as a distance between the end-situated heavy particles Req = R+ + R-, see Table IVTT|l . The 
lifetimes of these systems are very short. The ion {H ~ He — H)'^^ is highly unstable towards 
decays (HeH)'^^ + p, {He)~^ + 2p, but it does not decay to + a or H + a + p. In turn, 
the ion {He — H — He)^^ is highly unstable towards decays [HeHY'^ + a, He^ +p, but it 
does not decay to iJe"*" + a + p 01 H + 2a. 

Conclusion. We presented theoretical arguments about a possible existence of exotic ions 
{HeHf+ and Hef^ as well as some indications to the existence of the exotic ions Lir7 , 
[H — He — H)^^, {He — H — HeY^ . A striking conclusion can be drawn that the exotic 
molecular ion He^ is the most bound one-electron system among those made from protons 
and/or a— particles at strong magnetic field B > 1000 a.u. It is then followed by two 
other a-particle containing ions: the atomic ion He~^ and the hybrid ion [HeHY^ - are 
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characterized by slightly larger total energies (see Tables I, IV, VII). It is worth mentioning 
that iJ-atom - the only neutral one-electron system is the least bound one. Both ions He^'^ 
and {HeHY^ can exist also in the excited states. In Table lVni a comparison of total energies 
of different one-electron systems at the magnetic field B = 4.414 x 10^'^ G is given. In very 
strong magnetic fields 5 > 3 x 10^^ G the ions H, H^, H^^, He^, {HeHY^, He\^ are stable. 
In this domain of magnetic fields a striking relation between the binding energies of He^ 
and HI'^ exists: their ratio does not depend on magnetic field and is equal to two. It 
may allow to explain the observed ratio of the absorption features of the isolated neutron 
star 1E1207. 4-5209 0] by the helium-hydrogen content of its atmosphere under a surface 



magnetic field 4.4 x 10 G 
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